MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  bure.au  rn  STANDARDS  i%<  a 


OTIC  FILE  COPY  I  2  99  16 


Xi lo  /Ui 


Cornell  University 


LECTURES  ON  MATHEMATICAL  COMBUSTION 
Lecture  10:  Free-Boundary  Problems 
Technical  Report  No.  155 


J.D.  Buckmaster  &  G.S.S.  Ludford 


January  1983 


Theoretical  and 
Applied 

Mechanics 


DTlC.  ; 

JON  3011*3  ^ 


Thurston  Hall 
Ithaca ,  New  York 

88  06  30  010 


LECTURES  ON  MATHEMATICAL  COMBUSTION 
Lecture  10:  Free-Boundary  Problems 
Technical  Report  No.  155 

J.D.  Buckmaster  &  G.S.S.  Ludford 
January  1983 


U.S,  Army  Research  Office 
Research  Triangle  Park,  NC  27709 

Contract  No.  DAAG29-81-K-0127 

Cornell  University 
Ithaca,  NY  14853 


Approved  for  public  release;  distribution  unlimited. 


The  view,  opinions,  and/or  findings  contained  in  this  report 
are  those  of  the  authors  and  should  not  be  construed  as  an 
official  Department  of  the  Army  position,  policy  or  decision, 
unless  so  designated  by  other  authorized  documents. 


Contents 


Page 


1.  The  Hydrodynamic  Limit  1 

2.  The  Burke-Schumann  Limit  3 

3.  NEF  Tips  5 

4.  NEF  Wall -Quenching  8 

5.  Straining  NEFs  13 

6-  Shearing  NEFs  15 

Appendix.  The  Method  of  Lines  16 

References  19 

Captions  20 

Figures  1-9  21 


-10.1- 


Lecture  10 

FREE-BOUNDARY  PROBLEMS 

Throughout  these  lectures  we  have  ensured  that  the  reaction  terms 
vanish  everywhere  except  in  a  thin  (flame)  sheet,  whose  location  has 
to  be  found  as  part  of  the  solution.  So  far  this  free  boundary  has  been 
either  a  plane,  a  circular  cylinder,  a  sphere,  or  a  perturbation  of  one 
of  these;  we  now  consider  problems  with  more  complicated  free  boundaries  - 
There  are  four  different  ways  of  confining  the  reaction  to  a  sheet. 

(i)  Adopt  the  6-function  model  discussed  in  lecture  7- 

(ii)  Take  the  hydrodynamic  limit. 

(iii)  Take  the  Burke -Schumann  limit. 

(iv)  Use  activation-energy  asymptotics. 

In  this  lecture,  which  is  an  expanded  version  of  Buckmaster  (1982),  we 
shall  briefly  mention  examples  of  (ii)  and  (iii),  but  most  of  the 
discussion  will  deal  with  parabolic  problems  for  premixed  flames  that  arise 
from  (iv). 

1.  The  Hydrodynamic  limit 

Ca»o) 

The  Jump  conditions  ( 3 . 9); <j»  imply  that  the  flow  will  be  diffracted 
by  an  inclined  flame  S^r-ftte .  If  the  flame  speed  W  is  much  smaller  than 
the  speed  U  of  the  fresh  gas,  i.e. 

W  =  eU  with  e  «  1,  (l) 

a  uniform  flow 

y1  =  (u,o)  (2) 

is  turned  by  a  plane  flame  so  that  the  velocity  of  the  burnt  gas  is 


ye  =  (U,±eU(cr-l) )  +  o(e,e) 


(3) 
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where  the  sign  is  opposite  to  that  of  the  flame  slope  (figure  1).  The 
corresponding  jump  in  pressure  across  the  flame  is  0(e).  Such  a  deflection 
of  the  streamlines  is  a  fundamental  characteristics  of  tube  flames  (figure  2), 
which  are  often  slender  (i.e.  correspond  to  small  e). 

In  order  to  describe  the  shape  of  the  flame  when  the  prescribed  efflux 
from  the  tube  is  the  plane  flow  U(f(n),0),  we  look  for  i(nterior)  and 
e(xterior)  solutions 


u  = 


j  Uf (n)+eu^(ex,n)+. . .  je^v^ex  ,n)+. 

Uf (n)+eu^(ex,en)+. . . 


|ev*(ex,en)+.. 


(U) 


These  expansions  are  consistent  with  the  slender- flame  approximation, 
Euler's  equations,  and  the  jump  conditions.  Both  the  flame  shape  and 
the  flow  field  can  be  constructed  in  a  straightforward  manner  (Buckmaster 
&  Crowley  1982). 

In  fact,  a  differential  equation  for  the  shape 


n  =  ±F(ex) 


(5) 


i  2 

of  the  flame  can  be  deduced  immediately.  Since  v  is  0(e  ) ,  it  does 
not  contribute  to  the  0(e)  normal  velocity  ahead  of  the  flame,  so  that 
the  flame  speed  has  the  required  value  (l)  if 


1  +  f (F)F'  =  0, 


(6) 


a  result  that  is  also  true  for  axisymmetric  flow.  When 


f(n)  *  l-n2/a2  with  |n|  <_  a. 


(7) 


i.e.  the  flow  is  Poiseuille,  the  flame  shape  is  given  by 


F  -  F3/3a2  «  -ex- 


(8) 
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Only  for  slender  flames  can  the  shape  be  determined  directly  from  the 
given  flov  leaving  the  burner. 

2.  The  Burke-Schunann  Limit 

This  limit  arises  in  the  context  of  diffusion  flames  as  V  -*■  «.  As 
we  saw  in  section  8.  ,  the  fuel  is  then  absent  on  one  side  of  the  flame 
sheet  and  the  oxidant  on  the  other.  The  problem  considered  by  Burke  and 
Schumann  in  their  seminal  work  was  the  divided  flow  of  fuel  and  oxidant 
through  concentric  tubes  with  the  flame  attached  to  the  rim  of  the  inner 
tube  (figure  3).  If  the  oxidant  is  in  excess,  the  flame  terminates  at 
the  axis,  and  is  said  to  be  overventilated;  if  there  is  an  excess  of 
fuel,  termination  occurs  at  the  outer  tube,  an  underventilated  condition. 

Burke  and  Schumann  neglected  the  effect  of  the  flame  on  the  flow, 
i.e.  they  adopted  the  constant-density  approximation  over  50  years  ago. 

It  is  curious,  in  view  of  the  enormous  insight  provided  by  their  well- 
known  analysis,  that  to  this  day  one  finds,  at  scientific  forums,  objections 
to  the  approximation  when  it  is  used  ir.  activation-energy  asymptotics.  They 
also  neglected  longitudinal  diffusion,  a  step  that  can  be  legitimized 
by  setting 

X  =  x/U  (9) 

and  letting  U  *►  The  equations  to  be  considered  are  then 

3T/3x  =  32T/3y2,  3X/3X  =  K_132X/3y2,  3Y/3x  =  L_132Y/3y2  (10) 

under  the  jump  conditions 

6(l)  =  o(X)  =  5(Y)  =  6(3T/3y  +  2K_13X/3y)  =  6(3T/3y  +  2L_13Y/3y)  =  0.  (ll) 

The  plane  version  of  the  equations  used  by  Burke  and  Schumann  has  been 
written  because  we  shall  consider  a  simpler,  but  related,  problem  than 
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theirs.  General  Levis  numbers  have  also  been  introduced,  which  necessitates 
a  slight  modification  of  the  jump  conditions  derived  in  section  8. 

Consider  a  plate,  coincident  with  the  negative  x-axis,  separating 
parallel  and  equally  fast  flows  of  oxidant  and  fuel  ( figure  h) .  The 
initial  conditions  are 


r 

Y  =<  for  y  *  0  at  x  =  0; 


in  fact,  X  is  required  to  vanish  everywhere  below  the  flame  sheet  and 
everywhere  above.  The  solution  of  the  problem  (10-11)  can  be  written 
in  terms  of  the  similarity  variable 


(12) 


0! 


(13) 


as 


’x  L  erfc(^u/2)  ' 
‘L  erfc(XJu#/2) 

0 


Y  = 


r  _  erfc(I^ii>/2)  I 
'L  erfc(iJ0-/2)J 


for  u i  £ 


(1U) 


is  given  by 


2 

X_„L* s~K^/k  erfcHii»/2)  =  Yf^e"L“*Aerfc(^ V2)  - 


(15) 


Corresponding  formulas  for  T  could  also  be  written. 

"nen  the  reactants  are  supplied  in  stoichiometric  proportions  (i.e. 

X,  s Y^)  and  in  addition  K  =  l,  the  flame  sheet  coincides  with  the  positive 
v-axis,  i.e.  =  0.  A  decrease  (increase)  in  Xf  or  an  increase  (decrease) 
in  K  moves  nhe  flame  sheet  up  (down). 
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Our  simple  problem  has  yielded  to  analytical  treatment,  but  in  general 
the  Burke-Schumann  limit  involves  numerical  integration  of  a  type  that 
arises  in  the  free -boundary  problems  of  activation-energy  asymptotics 
to  be  treated  next . 

3.  KEF  Tins 

The  most  interesting  examples  of  free -boundary  problems  uncovered  by 
activation-energy  asymptotics  concern  ITEF's,  and  discussion  of  these  examples 
have  almost  invariably  adopted  the  constant-density  approximation.  The 
relevant  equations  and  jump  conditions  have  been  developed  in  section  4. 

The  starting  point  is  a  steady  plane  deflagration  in  a  uniform  flow 
with  speed  U  greater  than  1  (figure  5).  The  solution  is 

f  n  n 

T^.  +  ~IY  nen 

T  =  *  *  ,  h  ={  for  n  £  0,  (16) 

lHf  s  Tb  [  0 

formulas  that  have  already  been  used  in  connection  with  plane-flame  stability 
(section  5).  How  consider  the  effect  of  introducing  an  adiabatic, 
non-catalytic  wall  at  the  plane  y  =  0.  For  simplicity  we  shall  take 
i  =  0,  so  that  h  vanishes  identically  and  the  problem  reduces  to  one  for 
T  alone,  namely  the  free-boundary  problem 


U37/3x  =  V2T  for  n  <  0  and  y  >  0,  (17) 

7  -*■  T.  ,  3T/3n  -*■  Y.  as  n  -*■  0- ,  (l8) 

0  I 

3T/3y  =  0  at  y  =  0,  (19) 

T  -*•  7^  +  Yfen  as  x  -*•  -*  with  n  fixed.  (20) 


The  last  condition  comes  from  the  requirement  that,  far  from  the  wall,  the 
temperature  field  has  the  description  (l6a);  this  structure  gives  way  to  a 
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tvo-dinen s ior.al  combustion  field  as  the  wall  is  approached.  Here  we  have 
the  problem  of  a  plane-flame  tip  for  L  =  1,  the  wall  corresponding  to  the 
line  of  symmetry.  No  attempt  has  been  made  to  solve  it  in  this  form  because 
of  its  elliptic  nature. 

If  U  is  large,  the  formulation  becomes  parabolic,  giving  a  classical 
problem  of  Ste-for*  type.  Thus,  introducing  the  coordinate  (9)  and  letting. 
U  f  ®  converts  the  problem  ( 17-20)  into 


3T/3X  =  32T/3y2 

for  0  <  y  <  F(x), 

(21) 

Tv  3T/3y  -  Yf 

as  y  ->  F(x)-0, 

(22) 

3T/3y  =0  at 

n 

o 

<# 

(23) 

Tf  +  Yf  ey+x  as 

x  -*• 

(24) 

where 

y  *  F(x)  (25) 

is  the  free  boundary  and  the  origin  has  been  taken  at  the  intersection  of 
the  'undisturbed  flame  with  the  wall,  i.e. 

lin  f F ( x ) +X ]  *  0*  (26) 

x-*-00 

This  problem  is  identical  to  that  of  a  plane  flame  approaching  a  parallel 
adiabatic  wall,  with  x  playing  the  role  of  time. 

The  solution  of  the  Stephan  problem  (21-24),  and  similar  ones 
that  are  discussed  later,  must  be  obtained  numerically.  Meyer's  (1977) 
method  of  lines  is  narticularly  well  suited  to  this  task  and  is  outlined 
in  an  appendix.  It  was  used  to  obtain  figure  6,  which  reveals  a  significant 
increase  in  flame  speed  as  the  flame  sheet  approaches  the  axis,  i.e.  near 
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the  tip.  This  adjustment  on  the  diffusion-length  scale  of  a  TIE?  smooths 
out  the  sharp  tip  predicted  by  the  hydro  dynamic  analysis  in  section  1. 

A  similarity  solution  can  be  constructed  for  such  problems  in  the 
neighborhood  of  the  tip.  Setting 


T  =  +  (xt-x)JG(i)  with  u  =  (27) 

-here  x^  defines  the  location  of  the  tip,  leads  to  the  differential  equation 


G"  -  juSG'  +  jG  =  0 


and  the  boundary  conditions 


G'(0)  =  0,  G(ai#)  =  0,  G'Coj*)  =  Y  . 


Here  the  constant  <j#,  which  has  to  be  found,  determines  the  locally  paraboli 
flame  sheet.  The  solution  of  equation  (28)  satisfying  the  conditions  (29b, c) 


G  =  a  f  ex  dx/x  ; 


condition  (29a)  is  then  satisfied  ii 


*'  /?  o 

J  *  ex  dx  =  l,  i.e.  10^  =  1.85. 


The  similarity  solution  breaks  down  in  the  neighborhood  of  the  tip  where 
the  flame  slope  is  Q(l),  but  the  extent  of  this  neighborhood  is  given  by 

X.  -  X  =  0(u-2)  (2 


ar.d,  hence,  can  be  made  arbitrarily  small  by  taking  U  large  enough. 


-10.8- 


For  l  r  0,  i.e.  L  #1,  the  function  h  no  longer  vanishes  and  the 
problem  (21-2*0  must  be  augmented  with  equations  and  conditions  for  h, 
in  oarticular 


3h/3y  =  0  for  y  =  0.  (33) 

The  effects  of  Lewis  number  are  also  shown  in  figure  6.  Decreasing 
from  1  causes  the  tip  to  elongate  and  eventually  assume  a  bulbous  form. 

In  is  tempting  to  terminate  solutions  of  the  latter  type  at  the  point  Q 
in  figure  6,  for  they  then  closely  resemble  the  open  flame  tips  seen  in 
the  combustion  of  lean  hydrogen  or  rich  heavy-hydrocarbon  mixtures, 
mixtures  for  which  the  effective  lewis  number  is  significantly  less  than  1. 
3ut  the  mathematics  gives  no  clear-cut  reason  for  doing  this;  the  only 
suggestion  is  the  marked  decrease  in  temperature  along  the  flared  portion 
of  the  flame. 

The  negative  flame  speeds  associated  with  those  portions  of  the  flame 
with  positive  slope,  although  curious,  do  not  violate  the  physics.  There 
is  a  diffusive  flux  of  reactant  in  the  y-direction  towards  the  flame  sheet 
to  maintain  the  combustion. 

1.  '.7Z?  Wall-Onenchir.z 

When  the  wall  is  not  adiabatic,  the  problem  is  no  longer  relevant  to 
flame  tips.  For  a  cooled  wall  it  models,  in  a  rough  sense,  the  behavior 
of  a  flame  near  a  burner  rim  (figure  7).  It  corresponds  more  precisely 
to  the  propagation  of  a  plane  wave  towards  a  parallel  cooled  wall,  but 
our  discussion  will  be  couched  in  terms  of  burner  flames. 

A  flame  can  be  stabilized  on  a  Bunsen  burner  only  under  a  limited 
range  of  conditions.  Outside  this  range  either  "blow-off"  will  occur 
so  that  the  flame  becomes  detached  from  the  burner,  or  "flashback" 
will  take  place,  the  flame  traveling  to  the  base  of  the  burner  via 
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the  inner  wall  of  the  tube.  The  tendency  fcr  flashback  is  easy  to  under¬ 
stand:  a  prefixed  flame  will  travel  upstream  unless  the  gas  speed  is 
greater  than  the  flame  speed,  and  at  the  surface  of  the  burner  tube 
the  gas  speed  falls  to  zero.  Flashback  is,  therefore,  inevitable  unless 
there  is  some  mechanism  to  prevent  the  flame  from  reaching  the  surface. 

Heat  transfer  from  the  gas  to  the  tube,  quenching  the  reaction  that 
sustains  the  flame,  is  commonly  regarded  as  one  such  mechanism;  the 
free-boundary  problem  under  discussion  illustrates  this  mechanism. 

The  boundary  conditions  at  the  wall  are  taken  to  be 

3T/3y  =  k(T-Tf)/0,  3Y/3y  =0  at  y  =  0,  (3*0 

i.e. 

3T/3y  =  0,  3h/3y  =  k(T-Tf)  at  y  =  0  (35) 

in  the  TIE?  formulation.  At  the  same  tine,  the  problem  will  be  generalized 
slightly  by  considering  a  not  necessarily  uniform  flow 

Y  =  (Uf(y),0)  as  U  +  «.  (36) 

This  is  more  realistic  than  a  uniform  flow  when  the  no-slip  condition  at 
the  wall  is  satisfied  by  taking  f(0)  =  0.  We  shall  set  l  =  0  so  as  to 
filter  out  the  enthalpy  loss  or  gain  at  the  flame  due  to  unbalanced  diffusion 
of  temperature  and  reactant. 

The  problem  is  again  numerical,  but  certain  features  can  be  derived 
analytically.  So  long  as  the  flame  intersects  the  wall,  a  similarity 
solution  describes  its  behavior  there;  when  it  extends  to  x  =  ®>  there 
is  an  asymptotic  solution. 

Consider  first 


f(y)  =  1  * 


(37) 
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later  we  will  briefly  discuss  the  more  realistic  choice  (55).  We  shall  assume 
that 

T  \  as  x  '  (38) 

the  approach  being  exponential  in  and  that 


1 

x"JF(x)  0  as  x  -*■ 


(39/ 


where  F  is  the  free-boundary  function  (25).  Self-consistency  of  the 
resulting  asymptotic  description  will  be  checked  in  due  course. 

We  start  with  the  problem  for  h  (which  has  no  jumps  at  the  flame 
sheet)  and  write 


h 


+  h 


(2)  +  h(3) 


(U0) 


where  *n^,h^,h^  satisfy  the  heat  equation  (21)  individually,  the 
boundary  conditions 


3h(l)/3y  =  kYf,  3h(2)/3y  =  0,  3h(3)/3y  =  k(T-Tf-Yf)  at  y=0,  (Ul) 

and  vanish  exponentially  repaidly  as  y  -*•  The  condition  (38)  ensures 
(3) 

tr.at  h  is  exponentially  small  as  X"*“  and  so  need  not  be  considered 
further.  The  asymptotic  behaviors  of  the  two  remaining  functions  are 
determined  by  similarity  solutions  of  the  heat  equation.  Thus 

h^~  X5G(w)  with  u  =  y/xS ,  G  =  kYf [uerfc(u/2)-(2/*4f )e~“  (J»2) 

-  (2)  . 

ar.c.  h  is  asymptotically  a  sum  of  the  eigensolutions 

i  2 

X-5_4E(i)  with  S  =  e-^  ^  D2n(u//2);  (1*3) 


here  n  is  a  non-negative  integer  and  D  is  the  parabolic  cylinder 
function.  Since  all  the  eigensolutions  vanish  as  x  "*■  "»  we  need  not 
consider  h  ^  ^  further . 
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Only  h^  is  left  and  from  it  we  find 

i 

h/x  5  =  -2kYr,/V~  +  0(to)  for  a  snail .  (44) 

*1 

In  view  of  the  hypothesis  (39),  setting  a  =  x~??  ia  this  expansion  gives 


the  value  of  h  at  the  flame  sheet ,  so  that 

♦  ,  *  (2kYf/4T2)x*[l+o(l)]  as  k  +  ®.  (45) 

Consider  now  the  asymptotic  behavior  of 

T  =  Tb  +  YfT.  (46) 

The  appropriate  modification  of  the  problem  (21-214-)  is 

3T/3X  =  S^/Sy2  for  0  <  y  <  F(X),  (itT) 

T  v  0,  3T/3v  -  exp(-**/2)  asy+  F(X)-0,  (48) 

3T/3y  =  0  at  y  =  0;  (1»9) 

the  initial  conditions  (at  x  *  -°°)  are  omitted.  The  solution  is 
i  _  i  o  i  i 

F(x)  '  <*>#x\  T-  exp(-kYfX?/ At^)x'G(<3)  with  cD  =  y/X“,  (5d) 


where 

G"  +  a2G  =  0,  G'(3)  =  0,  G(u„)  =  0,  G’fZJ  =  1  with  a2=kYf/2/irT2.  (51) 

Tote  that  the  result  (50a)  is  consistent  with  the  assumption  (39)  under 
which  it  was  obtained.  The  constant  ,  which  has  to  be  bound,  devormines 
the  asymptotic  shape  (50a)  of  the  flame  sheet.  The  solution  of  the 
equation  (51a)  satisfying  the  conditions  (51b, d)  is 

G  =  -cosai/a  sin  ai#;  (52)- 
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the  condition  (51c)  is  then  satisfied  if 

aw*  *  (n+j)i7  (53) 

for  sons  integer  n.  A  further  req.uireT.ent  cor.es  from  the  maximum  principle 
for  the  heat  equation,  namely  T  <_  7^  everywhere .  The  only  acceptable  n  is 
zero  and  we  have 

5.  =  *5/k  (5U) 

A  similar  analysis  is  possible  when 

f(y)  =  y»  (55) 

a  more  realistic  choice  physically.  Now 

h(1)  -  x1/3g(w)  with  w  =  y/x1/3,  G  =  kYft3r(-l/3,u3/9)/r(-l/3),  (56) 

so  that 

h/X1/3  -  -32/3kYf/r(2/3)  +0(u)  for  u  small;  (57) 

and 

0*  =  (32/3kYf/r(2/3)7^)X1/3tl+o(l)]  as  x  +  “  (58) 

provided 

-l  /? 

X  /Jy(x)  ■+•  o  as  X  -*■  "•  (59) 

The  7-problem  has  the  solution 

7(x)  -5.X279,  ex?(-32/3kYfx1/3/2T(2/3)T^)x2/9G((u)  with  u  »  y/'x279,(6o) 

where 

G"+a213  «  0,  G'(0)  =  0,  G(u#)  =  0,  G'(u#)  =  1  with  a3=3~1/3kYf/2r(2/3)‘I^  .(6l) 
Finally,  the  maximum  principle  leads  us,  as  before,  to  the  conclusion  that 
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u5J_1 y3(2ai^^2/3) 
au^Jg/s  ( 2au)^2/3 ) 


u# 


r(2/3)T^r2 

2kYT, 


) 


1/3 


(62) 


where  r  is  the  smallest  zero  of  J  ,^(r)* 

These  results  suggest  a  simple  dichotomy:  for  sufficiently  small  heat 
loss,  the  flame  eventually  intersects  the  wall  (a  necessary  condition  for 
flashback)  •,  but  for  larger  values  of  k  the  flame  ultimately  moves 
away  from  the  wall,  in  particular  according  to  the  result  '(60a)  for 
the  flow  (55),  and  flashback  cannot  occur.  Figure  8  shows  Buckmaster's 
numerical  solutions ;  clearly  an  increase  in  k  tends  to  reduce  the  flame 
speed  near  the  wall.  For  moderate  values  of  k,  he  obtained  slightly 
bulbous  shapes  (cf.  flame  tips  for  L  <  l),  suggesting  that  this  would 
be  the  case  for  all  sufficiently  large  k.  ?Tunerical  re-examination  of 
question  has,  however,  provided  good  evidence  that  the  dichotomy  mentioned 
above  does  indeed  occur. 


5 .  Straining  :TEFs 

So  far  we  have  dealt  only  with  parallel  "flows ,  but  parabolic  free¬ 
boundary  problems  can  also  be  formulated  for  the  more  general  velocity 
fields 

v  =  Ua(x/U,y/U)  with  U  »  1.  (63) 

These  are  fast  flows  with  Yy  =  0(l);  more  precisely,  each  is  effectively 
made  up  of  an  uniform  flow  and  a  simple  strain,  both  of  which  vary  slowly. 
The  angle  between  the  flame  and  the  streamlines  is  small,  so  that  it  stays 
close  to  a  single  streamline.  If  x  is  measured  along  the  streamline  and 
y  perpendicular  to  it,  the  velocity  in  the  neighborhood  of  the  stream¬ 
line  is  approximated  by 


Y  =  (UqQ,-q^y) 


(64) 
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v'r.ere  qQ(x)  is  speed  on  y  =  0.  In  the  limit  U  -*■  «•  the  NEF  equations 
(1.  )  become 

(qo3/3X-q^y3/3y)(T,h)  =  32(T.h+lT)/3y2,  (65) 

and  these  are  amenable  to  numerical  treatment  for  any  choice  of  a  .  Modest 

o 

variations  in  l  do  not  have  a  substantial  qualitative  effect,  so  that  ve 

* 

shall  set  l  =  0. 

Two  examples  will  be  considered,  namely 

qQ(x)  -  *35  sinx  for  0  <  x  <  it  ,  .35(x-2x3/3-x5/5)  for  (66) 

Each  of  these  has  a  stagnation  point  as  x  0,  where  the  analysis  leading 

to  equations  (65)  breaks  down;  nevertheless  the  corresponding  stagnation- 

point  solution  constructed  in  section  k.  provides  the  correct  initial 

conditions.  Figure  9  shows,  for  each  example,  a  correlation  between 

the  flame  sreed  and  the  strain  rate  a',  a  measure  of  the  stretch 

"o 

experienced  by  an  element  of  the  flame:  as  the  stretch  increases  (decreases) 
the  flame  speed  decreases  (increases). 

The  conclusion  is  not  necessarily  valid  in  other  circumstances, 
however,  as  can  be  seen  from  analytical  results  obtained  by  Buckmaster 
(1532)  for  a  general,  small  but  rapid  velocity  change 

qQ  =  1  +  eQ(x)  with  0  <  e  <<  1,  x  =  x/e2  (67) 


on  the  scale  of  x*  The  wave  speed  has  the  explicit  representation 


dy. 


(68) 


a  formula  that,  surprisingly  enough,  is  valid  for  all  values  of  t.  There 
is  no  strong  connection  between  the  behavior  of  W  and  changes  in  the  stretch 
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eQ* i  in  fact,  a  change  in  the  sign  of  the  stretch  has  more  influence 
on  V. 

The  result  (63)  shovs  that  Q  can  be  chosen  (vith  a  degree  of 
arbitrariness)  to  mahe  W  vanish  at  any  given  point.  The  possibility 
of  decreasing  the  flane  speed  to  zero  at  a  point  by  adjusting  the  velocity 
field  has  not  been  demonstrated  before. 


6.  Shearing  !T5Fs 

Me  conclude  our  discussion  of  free-bourdary  problems  by  briefly 
describing  the  effect  of  shear  on  a  premised  flame.  The  vail  is 
removed  and  a  linear  shear  flov  is  inserted  for  y  <  0,  the  flow  being 
assumed  uniform  for  y  >  0.  More  precisely, 

fl 


f(y)  - 


l-uy 


for  y  £  0 


(69) 


in  the  expression  (36). 

'Numerical  results  have  been  obtained,  in  particular  for  u>  =  5.  Here 
ve  shall  be  content  to  describe  the  solution  in  the  limit  of  very  strong 
shear,  i.e.  <»>-*■<*>.  Then  the  equations  in  the  lower  half-plane  simplify 


to 

3T/3X  =  3h/3X  =  0, 


(TO) 


giving  the  solution 

T  =  Tf,  h  =  0;  (71) 

these  become  boundary  conditions  (  at  y  =  0)  on  the  solution  in  tie  upper 
half-plane . 

The  asymptotic  behavior  as  x  ®  is  particularly  simple,  depending 
only  on  y : 
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T  = 


W /F(") 


h 


T. 


-vc^f  ?(») 


-Vi. 


for  y  <  F(«), 


vr.ere 


F(») 


l 

e  . 


(72) 


(73) 


numerical  integration  shows  how  the  renote  combustion  field  (16)  is 
transformed  into  the  asymptotic  field  (72).  The  flame  speed  vanishes 
as  x  ■*  ”  and  is  negative  or  positive  in  the  neighborhood  of  infinity 
accordingly  as  i,  is  positive  or  negative. 


Appendix.  The  Method  oTf  Lines 

Consider t’ e  parabolic  free-boundary  problem  (21 -2U).  Let 

Tn(y)  =  T(xn.y),  K  =  F(x _).  (7M 

ar.d  approximate  the  partial  differential  equation  (21)  by  the  ordinary 
differential  equation 

T"  =  T  /cx  -  ?  ./Ax  with  Ax  *  X  -X_  (75) 

n  n  n— j.  n  n— i. 

Having  determined  7  1  at  the  previous  step,  we  must  integrate  this 

second-order  equation  for  T  along  the  line  x  =  X  subject  to  three 

n  n 

boundary  conditions  (22,23).  The  integration  will,  therefore,  determine 

non  only  T  but  also  ?  . 

J  n  n 

Let  A(y)  be  the  solution  of  equation  (75)  with  A(0)  =  1,  A'(0)  =  0 
and  let  3(y)  be  the  solution  with  A(0)  =  A'(3)  =  1.  Then 

T^  =  (/-3)b  +  3  (76) 

is  the  solution  satisfying 

T  (0)  =  o,  T’(0)  =  0, 
n  n 


(77) 
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ar.d  from  it  we  may  calculate 


T*  =  (A'-B')o  +  3'.  (78) 

n 

Plimir.ation  of  b  from  the  equations  (75), (?S)  gives  a  relation  of  the 
forn 

T*  =  ?T  +  Q,  (79)" 

r.  11 

where  ?(y)  and  Q{y)  could  be  written  in  terms  of  A  and  B.  The  initial 
conditions  (77)  show  that 

P(0)  =  Q(0)  =  0,  (80) 

since  ?,Q  (like  A,B)  are  independent  of  the  parameter  b.  The  key  step 

is  to  find  differential  equations  for  P,Q. 

If  T'  and  T"  are  eliminated  between  the  two  equations  (75), 
n  n 

(79)  we  are  left  with 

T  (P'+P2-1/Ax)  +  (Q'+PQ+T  ./Ax)  =  0.  (8l) 

n  n— x 

Since  only  T  denerds  on  b,  the  parentheses  must  separately  vanish, 
n 

giving  a  pair  of  first-order  differential  equations  for  P,Q.  Numerical 
integration,  under  the  initial  conditions  (80),  then  determines  these 
functions . 

The  nosition  y  =  ?  of  the  free  boundary  can  now  be  determined 
n 

as  a  root  of  the  equation 

Yf  =  +  Q(Fn>»  (82) 


obtained  by  substituting  the  flame-sheet  conditions  (22)  in  the  relation 
(79)-  With  F  found  in  this  way,  the  relation 


(considered  as  a  first-order  differential  equation  for  Tn)  can  be  integrated 

towards y  =  0  with  the  initial  condition  T  (?  )  **  T.  ,  to  determine  T  (y). 

n  n  b  n 

The  method  is  quick,  efficient,  and  requires  no  iterations. 
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(i)  Adopt  the  6-function  model  discussed  in  lecture  7. 

(ii)  Take  the  hydrodynamic  limit. 

(iii)  Take  the  Burke- Schumann  limit. 

(iv)  Use  activation-energy  asymptotics. 

In  this  lecture,  which  is  an  expanded  version  of  Buckmaster  (1982),  we  shall 
briefly  mention  examples  of  (ii)  and  (iii),  but  most  of  the  discussion  will  deal 
with  parabolic  problems  for  premixed  flames  that  arise  from  (iv). 
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